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Abstract. We study density of states and conductivity of the doped double-exchange system, treating 
interaction of charge carriers both with the localized spins and with the impurities in the coherent poten- 
tial approximation. It is shown that under appropriate conditions there is a gap between the conduction 
band and the impurity band in paramagnetic phase, while the density of states is gapless in ferromag- 
netic phase. This can explain metal-insulator transition frequently observed in manganites and magnetic 
semiconductors. Activated conductivity in the insulator phase is numerically calculated. 

PACS. 75.50.Pp Magnetic semiconductors - 75.30.Vn Colossal magnetoresistance - 72.10.-d Theory of 
electronic transport; scattering mechanisms 

1 Introduction The present paper is devoted to the consideration of 

single-particle states and conductivity in impure DE sys- 
tem. Interaction of charge carriers both with the local- 
The recent rediscovery of colossal magnetoresistance (CMR) ized spins and with the impurities is strong , so it is defi . 

m doped Mn oxides with perovskite structure Ri_ I D I Mn0 3 nitely not enough to limit ourselves with the finite number 
(R is a rare-earth metal and D is a divalent metal, typ- of terms of per turbation expansion. A simple but physi- 
ically Ba, Sr or Ca) W has generated substantial inter- cally mea ningful approximation, allowing to sum up in- 
i est in these materials g. The doping of parent material finite number of perturbation expansion terms is the co- 
in RMnOa by a divalent metal is the source of the holes re- herent potential approximation (CPA). Initially CPA was 
sponsible for the transport properties of these materials. proposed to treat potential disorder @, but soon after it's 
In addition, each divalent atom introduced, is the center appea rance the generalization to random spin system was 
of an impurity potential. Many papers analyzed the mflu- dev eloped @. The CPA was also used to describe diluted 
ence of strong magnetic disorder, inherent in the the CMR magnetic semiconductors S 

materials at finite temperature, upon the single-particle In the present paper wc for thc first time treat Qn 
states and transport properties. However, the interplay equ&1 footmg the interactions of electrons with the core 
between the magnetic disorder and the dopmg-mduced Mn spins and with the doping impur ities using the matrix 
disorder was studied less. The impurity potential plays generalization of the CPA. The concurrent action of po- 
double role. First, the potential fluctuations determine tential disorder and tem perature dependent spin disorder 
the transport at temperatures well below the ferromag- leadg to a number of interesting phenomena, in particular 
net (FM) - paramagnet (PM) transition point T c . Second, to the possibility of the opening of the gap at the Fermi 
strong potential may pm the Fermi level either m the con- levd with the increase of temperature and, hence, to MIT 
duction band tail (in the Anderson model of disorder ||3Jj ) , transition 
or in the emerging impurity band. The analysis of exper- 
imental data reveals strong relevance of the latter effect 
to metal-insulator transition (MIT) near T c both in mag- 
netic semiconductors {§ and manganites [§. However, to 2 Hamiltonian and Theoretical Formulation 

the best of our knowledge the impurity-band scenario in 

the double-exchange (DE) model was not discussed yet. We consider the DE model with the inclusion of the single- 

site impurity potential. In addition, as it is widely ac- 

a e-mail: marka@ee.bgu.ac.il cepted, we apply the quasiclassical adiabatic approxima- 

b e-mail: kogan@quantum.ph.biu.ac.il tion and consider each Mn spin as a static vector with 
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a fixed length S (Si = ^n;, where is a randomly ori- 
ented unit vector). The Hamiltonian of the model in site 
representation is 



ti-j + Sij (tj — Jlii ■ a) = Hkin + Vi m p + V s d, (1) 



where ti-j is the electron hopping, is the random on-site 
energy, J is the effective exchange coupling between a Mn 
core spin and a conduction electron and a is the vector of 
the Pauli matrices. The hat above the operator reminds 
that in one-particle representation it is a 2 x 2 matrix in 
the spin space (we discard the hat when the operator is a 
scalar matrix in the spin space). 
We present Hamiltonian as 

H = + E + V imp + V sd ~E = H Q + V (2) 

(the site independent self-energy E(E) is to be determined 
later), and construct a perturbation theory with respect 
to random potential V — Vi mp + V s d — E. To do this let 
us introduce the T-matrix as the solution of the equation 



where 



T = V + VG T, 
G» 1 



E-H Q 

For the exact Green function we get 
G = Gq -+- GqTGo 



(3) 
(4) 

(5) 



The coherent potential approximation (CPA) is expressed 
by the equation 

'&) = G . (6) 



This equation can also be presented as 

'fi) = : II. 



where Tj is the solution of the equation 



and 



f i = V i + V i g(E-E)f i , 



g(E) = (G (E)) U = J ^Mcfe, 



(7) 



(8) 



(9) 



where Nq(e) is the bare density of states. The averaging in 
Eqs. (^,0) should be performed both with respect to ran- 
dom orientations of core spins and with respect to random 
on-site energies. We obtained, in fact, the algebraic equa- 
tion for the 2x2 matrix E 



1 - V i9 (E - E) 



Vi 



0. 



(10) 



This equation takes into account scattering both due to 
randomness of the core spins, and due to the impurities. 
If the impurity potential is negligible (V = 0) this equa- 
tion coincides with the Eq.(20) of Ref. 0] obtained in the 



dynamical mean field approximation (and also with those 
obtained for the Falikov-Kimball model p||lT|| ). 

In the reference frame where the z axis is directed 
along the magnetization, E is diagonal, and Eq.([l0|) re- 
duces to the system of two equations for its diagonal ma- 
trix elements E a (E) (a =f, |). The equations acquire es- 
pecially simple form at two extreme particular cases, which 
we will analyze: 

(a) T — 0. The magnetic state is coherent FM with n| = 1, 
and Eq. (|Io|) takes the form 



e t =p J - S^ tl 



l-{eiTJ-Eu)9(E-E u 



= 0. 



(11) 



(b) T >T C and zero magnetic field. The magnetic state is 
isotropic PM with (rii) = 0, which leads to E± = = E, 
and Eq. (|Io|) takes the form. 



J — E 



l-{e i + J-E)g{E-E) 

Ei - J - E 



I- {a- J- E) g{E - E) 



= 0. 



(12) 



We will solve the equations Eqs.(jll]) and (|T^) in the 
strong Hund coupling limit (J — > oojT In this limit we 
obtain two decoupled spin sub-bands. The equation for 
the upper sub-band, after shifting the energy by — J, for 
both cases (a) and (b) can be written down in unified form 



1 



l-(e l -E)g{E-E) 



(13) 



where a — 1 for T — 0, a — 2 for T > T c . In the model 
of substitutional disorder (ej = with probability x, and 
e,; = V with probability 1 — x), Eq. (|l3|) takes the form 



1 - x 



l + Eg{E-E) 1+ (E-V)g(E-E) 



(14) 



3 The CPA equations for semi-circular bare 
density of states 



We consider semi-circular (SC) bare DOS 

4 



No(e) = —Ml 



(15) 



at \e\ < W/2 and N (e) = otherwise, for which 



— -,/(— I -1 

w V V w 



(16) 



Let us introduce the following normalized quantities 



E 
W 



to 



E 



v 
If 



(17) 
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After simple algebra we obtain from Eq. (14) the cubic 
equation with respect to 



7 = Wg (E - S) = 8 Co- 
in the form 

7 3 + 16 (v - 2lu) 7 2 + 16 



A -y/(ui- A) 2 -1/4 



,(18) 



16uj (v — uj) 

v 



7 

-256- + 256 (1 - x) - = 0. (19) 
a a 

The number of electrons per cite n is given by 

f(E)N(E)dE, (20) 



(21) 



where f{E) is the Fermi distribution function, and 
N(E) = — Im 7 



is the actual density of states. To define the position of fi, 
the Fermi level, we must impose the relation between n 
and x; the simplest assumption appropriate for mangan- 
ites is the equation n — 1 — x. 



4 Conductivity in CPA 

For a disordered one-electron system the static conductiv- 
ity is given by 



is the one-particle spectral weight function. On account 
of the locality of T-matrix the second term in the trace is 
equal to 

^ Ss'vkaVk'aGa (e k , E s ) G„> (e k ' , E„) 

s,s — ±; k,k' ,cr,<r' 

xG a i (eic',E s t) G a (ek, E s i) 
x {T aa ,{V- \s!,E a )T a > a {\s! -k,E s ,)). (27) 

Since in CPA (T aa ,(k - k', E s )T al(T {k' - k, E s ,)) does not 
depend on k and k' and v_k Q = —viaa the above expres- 
sion is identically zero |l^,[ll|. Thus, finally 

_ a e 2 irh f f ( Of 

■v 2 a {e)N (e)J2[Aa(e,E)] 2 dEde, (28) 
where v is the unit cell volume and by definition 

,S(e)JV ' (e) = W?(l) ,J(e - ek) " (29) 

Let us assume nearest-neighbor tight binding spectrum 
on simple d-hypercubic lattice (v — a d ) 



£k = —t cos ak a , 



a=l 



Tr 



e nh 
V 



v a S [E — H) v a 5 [E — H 



df_ 

dE 



v 2 a (e)No(e) 
In SC DOS model 

vl (e)N (e) 



eN (e) de. 



(30) 



4 a 2 



ttW dh 2 



z\ 1 



W/2 



dE, 



(22) 



1 W 



where V is the volume and v a is a Cartesian component of 
the velocity operator. To obtain the conductivity in CPA 
let us express operator delta-function as follows 



3 7r dh 2 V W 2 
Substituting this result into Eq.(J2q) we obtain 



4s 



2z 
W 

2\ 3/2 



dz 



1 

2vri 



Using Eq. (§) and Eq. (|) in Bloch representation 

<5k,k"5<T,<r' 



krr 



Go(^) 



kV 



£ - e k - S a (E) ' 



(23) 



(24) 



with 



00 



^0 



a/ 

<9£ 



we get 



Tr [v a S (E - i?) v a 5 (E-H 

= J2< (ek,E)] 2 + o((rr\), 



2Trda d - 2 h 

being the Mott minimal metallic conductivity, and 
3^ J_ w/2 \ VF 2 



(32) 
(33) 



where 



k.cr 



1 



A a (e, E) = — 



ImZv (E) 



Tr [E - e - ReE a (E)] 2 + [ImZV {E)f 



(25) 



(26) 



E-e-S a {E) 
For the strong Hund coupling we obtain 



de. 



(34) 



A(E) 



3tt 



1-a; 



2\3/2 





1 




|lm 


} 




a; — z 





dx, (35) 



where z = 2 (w — A). Eqs. (|32j), ( |35| ) give the conductivity 
in the framework of bare SC DOS model for arbitrary hole 
concentration x and impurity potential strength V. 
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Fig. 1. DOS at T = (full line) and high-temperature (dashed 
line) in the units of W' 1 for x = 0.2 and V/W = 0.4 



5 Influence of the impurity potential 

First, consider density of states. It is known, that within 
CPA for every x there exists critical value of potential-to- 
bandwidth ratio v c (x) such that at v > v c (x) the separate 
impurity band splits off the conduction band (that is a 
gap opens in N(E)). In our approach we get two differ- 
ent curves v c (x,T = 0) and v c (x,T > T c ), which present 
boundaries of metal-insulator and metal-semiconductor 
'phase diagrams', respectively, in (v, x) plane. Due to ef- 
fect of magnetic disorder it appears that v c (x,T — 0) > 
v c (x,T > T c ). 

For a typical concentration x — 0.2 v c (0.2,T = 0) w 
0.49 and u c (0.2,T > T c ) « 0.35. So if we choose v = 0.4 
both N(E) and A(E) must be gapless at T = but do 
have a gap at T > T c . Numerical calculations of the DOS 
performed at T = and T > T c for the above x and 
v clearly demonstrate FM-PM transition induced band 
splitting (Fig.l). 

Now address the question of conductivity. Consider 
first the position of /j, and conductivity at T = 0. We get 
from Eq.@ /i(T = 0) = 0.3662VF. Note that fi(T = 0) 
lies on the neck connecting conduction band and impurity 
states derived parts of the band. As a result, the residual 
conductivity (Eq.Q) p-\T = 0) = 0.6163o- is less then 
the Mott limit. 

At T > T c DOS and A(E) have the same gap A = 
0.031W (see Figs.1,2), so (i(T > T c ) must lie in the gap. 
Thus, the model describes a bad metal at T = 0, and a 
semiconductor at T > T c . The transition between two 
types of conduction (FM-PM transition induced MIT) 
should occur at some temperature below T c , Such a pic- 
ture agrees with the recent photoemission experiments 
showing drastic decrease of DOS at the Fermi level jjdj 
as temperature increases towards T c . 

It is checked numerically that DOS displays square- 
root like behavior near the top of the conduction band 
E n 



Fig. 2. Energy-dependent conductivity at T = (full line) and 
high-temperature (dashed line) in the units of <ro for x = 0.2 
and V/W = 0.4 



and the bottom of the impurity band E$ 
N(E)^n iy /E-Ei. 



(37) 



Unlike DOS A(E) behaves linearly near the band edges 

A{E) w W~ X \ C {E C - E), for E < E c ; 
A(E) w W~ 1 \ i (E - Ei), iovE>E l . (38) 

The assumption T < A allows us to explicitly obtain 
fi(T > T c ). Calculating integrals in Eq.pj) with expo- 
nential accuracy we obtain 



1 



E c + Ei+Tln- 



(39) 



The integral in Eq. (j32|) , calculated with the same accu- 
racy, leads to activation law for conductivity with linear 
temperature pre-exponent 



BT 



Ea 
T 



(40) 




N(E) mn c y/E c -E, 



(36) 



where Ea = A/2 « O.OlSly and B is the following nu- 
merical constant 



(41) 



for the parameters considered. 

Low values of conductivity obtained for the case of spin 
disorder are an indication of the possibility of Anderson 
localization [[|jlj,[l6|] , which CPA is incapable of accessing. 
But the present results complement and support the local- 
ization based approach. In fact, the results of Ref. B were 
obtained under the assumption of the Fermi level pinning, 
which is now explained as being due to strong electron- 
impurities interaction (and the impurity band formation). 

In another aspect, the model considered may also ex- 
plain low-temperature MIT observed in initially metallic 
manganites Ri^^D^MnOj upon substitution of R by iso- 
valent atoms (e.g. La by Y ]17||). One may speculate that 
the substitution forms a deep impurity band which can 
capture holes in Ri-^D^MnOa. 
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6 Conclusion 

To conclude, we derived CPA equations for the one-electron 
Green function and conductivity of DE system contain- 
ing impurities. The equations were solved for the SC bare 
DOS and substitutional disorder model. It was shown that 
if the electron-impurity interaction is strong enough, there 
is a gap between the conduction band and the impurity 
band in PM phase, the density of states being gapless 
in FM phase. Under appropriate doping conditions the 
chemical potential is pinned inside the gap. This can ex- 
plain metal-insulator transition observed in manganites 
and magnetic semiconductors. 
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